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ABSTRACT 



We present Spitzer 8 /xm transit observations of the extrasolar planet HD 
149026b. At this wavelength, transit light curves are weakly affected by stel- 
lar limb- darkening, allowing for a simpler and more accurate determination 
of planetary parameters. We measure a planet-star radius ratio of Rp/Ri, = 
0.05158 ± 0.00077, and in combination with ground-based data and independent 
constraints on the stellar mass and radius, we derive an orbital inclination of 
i = 85?4 and a planet radius of Rp = 0.755 ±0.040 Rj. These measurements 
further support models in which the planet is greatly enriched in heavy elements. 



Subject headings: stars: planetary systems — techniques: photometric 



1. Introduction 



Much attent ion has been lavished on the transiting extrasolar planet HD 149026b 
(ISato et al.l l2005l ) due to its potential to directly test models of planet formation. The 
planet's small observed radius for its mass imply t hat an extraord i nary fraction of its 



mass (roughly 2/3) is in the form of heavy elements (ISato et al.ll2005l : iFortney et al.l 12006 



Burrows et al.l 120071 ). The discovery of a metal- laden planet orbiting a very metal-rich host 
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star ([Fe/H] = 0.36: ISato et al.ll2005l ) strongly suggests that core-accretion (e.g., Pollack et 
al. 1996) plays a role in forming giant planets. If, however, most of the heavy elements 
reside in the planet's core, then HD 149026b would possess a core mass much gr eater than 
the expected critical core mass of 10 — 2OM0 (jMizundll98d : iPoUack et al.lll996l ). and thus 
nonetheless present a challenge to standard core-accretion theory. 

T he planet is noteworthy in another respect. Observations by the Spitzer Space Tele- 
scope ( iHarrington et al.l 120071 ) have shown the planet to have a day-side 8 /im brightness 
temperature well in excess of its predicted blackbody temperature, when it is assumed that 
all incident radiation is absorbed and sub sequently re-emitted uniformly across the entire 
surface of the planet. iFortney et al.l (120081 ) posit that highly irradiated planets such as HD 
149026b, which they term "pM" class planets, will generally show bright day-sides and large 
day/night temperature contrasts. They argue that the incident stellar flux is prominently 
absorbed by gaseous TiO and VO high in the atmospheres of pM plan ets where the radia- 



tive timescale is much shorter than the advective timescale (see also iHubeny et al.l 12003 



Burrows et al.ll2008l ). This is in contrast to less irradiated "pL" planets where Ti and V 
are expected to largely condense out of the atmosphere, permitting the stellar flux to be 
absorbed deeper in the atmosphere where the two timescales are comparable. Hence, it is 
only for the pL class that a heated parcel of gas is able to be advected to the night side prior 
to cooling, resulting in similar day/night temperatures. 

HD 149026b is thus a valuable case study for modelers of planetary atmospheres, struc- 
ture, and formation. Unfortunately, the system is observationally challenging: the transit 
depth (3 mmag in V) is a factor of two shallower than any other presently known transiting 
planet, and more importantly, there are few adequate comparison stars nearby on the sky. 
The result is that the pres ent fractional unc ertainty in the key observable parameter, the 
planetary radius Rp, is 7% (IWinn et al.ll2008l ). This uncertainty is one of the largest among 
the ensemble of transiting planets. The state of uncertainty is unfortunate given that Rp 
is the essential constraint on models of the planet's interior structure. Fortunately there is 
further scope for improvement through high-precision photometry. 

This study is inspired by the potential of infrared photometry with the Spitzer Space 
Telescope to reduce the uncertainty in Rp. While ground-based photometry suffers from 
significant levels of systematic noise when there are few good comparison stars, Spitzer has 
demonstrated 0.1 mmag photometry without any comparison stars (e.g. Knutson et al. 
2007). Additionally, because of the near absence of stellar-limb darkening in the infrared, 
transit light curve modeling is simplified and g ives results l argely independent of assumption s 
abo ut hmb-darke r iing c oefficients. Previously, ISato et al.l (120051 ) , ICharbonneau et al.l (120061 ) , 
and I Winn et al.l (120081 ) (hereafter W08) have presented ground-based photometry of HD 
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149026. In this paper, we report Spitzer 8 fim observations of the transit of HD 149026b, 
and combine this with the previously pubhshed data in order to derive precise constraints 
on the properties of HD 149026b. In §2 we describe the observations and data reduction 
and in §3 we describe our analysis of the Spitzer light curve. In §4 we estimate the physical 
parameters of the HD 149026 system. We conclude with a discussion of the implications of 
our revised estimate of the planet radius for models of the interior structure of HD 149026b. 



2. Observations and Reduction 

We observed the transit of HD 149026 on UT 2007 August 14, using the 8 micron 
channel of the IRA C instrument (Fazio et al 2004) aboard the Spitzer Space Telescope 



( IWerner et al.l l2004j ) . The system was observed at a 0.4 s cadence using IRAC's 32 by 32 
pixel sub-array mode, in which frames of 64 images are taken in rapid succession. Over the 
course of our observations, we obtained 1047 such frames, resulting in 67,008 total images. 
Our observational strategy matches that of recent Spitzer observations of HD 189733 and 
GJ 436 (e.g., Knutson et al. 2007, Deming et al. 2007, Gillon et al. 2007); the telescope 
positioning was held fixed to avoid time loss during telescope movements and to minimize 
errors from an imperfect flat-field correction. In the IRAC 8 micron channel, there is a 
well-known rise in detector sensitivity during observational sequences (see e.g., Harrington 
et al. 2007, Knutson et al. 2007), which is steepest at the beginning of observations and 
asymptotes within several hours for highly illuminated (> 250 MJy Sr~^) pixels. We padded 
the beginning of our observational window so that the transit would begin nearly 3 hours 
into observations, thus avoiding the steepest part of this "ramp." 

In each image, we assessed the background flux by taking the median pixel value from 
the corner regions of each 32 by 32 image. We performed aperture photometry, settling on 
a 3.5 pixel aperture radius, for which the rms of the time series is minimized. From the 
time stamp reported for each frame of 64 0.4 s exposures, we calculated the JD of the center 
of integration for each image. We applied the heliocentric correction to the JD using the 
position of Spitzer obtained from the JPL Horizons Ephemeris System. In each series of 64 
images there is a well-known effect, with the first 5-10 and 58*^^ images showing anomalously 
low star fluxes and background levels (see e.g. Harrington 2007). Background subtraction 
generally corrects for this effect, but we elected to drop the I'** image from each series of 
64, because the background levels in these images exhibit more dispersion than in the other 
images. We trimmed the first 45 minutes of data, when the ramp is steepest. We flagged 
images when the star centroid, calculated with a flux-weighted average, was 4a away from 
the median centroid position. Such centroid deviants were generally caused by cosmic 
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rays or other contamination in the photometric aperture. We further flagged images when 
the flux measurement was 4a from a smoothed (binned) hght curve, or the background level 
was 4(7 from a binned time series of the background. We flagged 317 images (0.5 % of the 
total) according to the last three criteria. 



3. Spitzer Light Curve Analysis 

One major benefit of observing transits at 8 fim is that stellar limb-darkening has a 
small effect on the shape of the tr ansit light curve. T o determine its extent, we consulted a 
theoretical limb-darkening model JKuruczlll979lll994)^ra T^s_ = 62507^. log^ = 4.5, [Fe/H] 



0.3 star at A = 8 /xm. We fit this models] to the Iciaret (2000) four parameter nonlinear 



limb-darkening law (see also Beaulieu et al. 2008 for a similar handling of limb-darkening). 
Though the limb-darkening is indeed modest, we find that incorporating it in our light 
curve modeling (described below) leads to non-negligible changes in the best-fit parameters 
and a reduction in the best-fit by more than 1. We modeled the light curve using the 
"small-planet" transit routine of Mandel & Agol (2002) . The small planet approximation 
is not usually suitable for analyzing high quality transit data, especially for systems with 
large planet-star radius ratios {Rp/ R^, > 0.1), but here we find the approximation leads to 
insignificant changes in the best-fit parameters (due to the very small planet-star radius ratio 
of HD 149026). We assumed a circular orbit, which is expected from tidal dissipation and 
supported by current radial velocity data (e.g., Sato et al. 2005). We parametrized the light 
curve with 4 geometric parameters that are independent of prior assumptions on the stellar 
properties: the planet-star radius ratio Rp/R^, the stellar radius to orbital radius ratio -R^/a, 
the inclination i, and the time of mid-transit T^. To correct for the ramp and other possible 
detector effects, we adopted a correction factor / = (cq + Ci log(t — to) + C2log^(t — to)), 
where to was fixed to a time a few minutes before the first observations. Note that in all of 
our modeling below, we fit for the detector correction coefficients simultaneously with the 
transit-related parameters, allowing us to take into account how changes in the correction 
coefficients may impact the transit parameters. 

We performed a least-squares fit to our unbinned data over the 7 parameter space 
{Rp lR-y, Ri,la,i, Tr., Co, ci, C2), using an IDL implementation of the amoeba algorithm (e.g. 



see 



Press et al.l Il992l ). The data, corrected for the ramp and binned 100:1 are shown in 
Figure [H together with this best-fitting solution. To understand the level of photometric 
noise and its properties, we examined the residuals from this best fit. We determined a 
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normalized rms residual of 8.3 x 10~^, only 15% greater than the expected photon-noise. We 
found that the level of photometric noise was constant over the duration of the observations, 
and furthermore that the noise was essentially "white." In the left panel of figure [2] we show 
that the scatter in binned residuals decreases with bin size as A^~^/^ for bins of up to 1000 
images. In the right panel of Fig. 2, we display a power spectrum estimate for the time- 
series of residuals. To compute this, we first binned the residuals for each frame of 64 images. 
This step creates an evenly spaced time-series of 940 residuals because IRAC, in sub-array 
mode, takes exposures in sets of 64 images (once every 25.6 seconds). This binning also 
avoids having to interpolate over gaps caused by flagged images. Though binning removes 
the highest frequency information from the spectrum, we are less concerned with noise power 
on the affected time scales, which are shorter than the other timescales relevant in a transit 
light curve (e.g., the ingress/egress duration). We estimated the power spectrum via the 
(modulus squared) discrete Fourier transform. To reduce the variance at each frequency, 
we smoothed with a 7 point "Daniell", or moving average, filter. We compared our power 
spectrum estimate with that expected of white noise by simulating 10^ time-series with 
identical, independent Gaussian deviates of the same time sampling and standard deviation 
as the Spitzer residuals. We note that 5% of the simulated power spectra show peak values 
as high as the peak value in Spitzer power spectrum, while two peaks in the Spitzer power 
spectrum exceed the median peak value of the simulated spectra (the dashed line in Fig. 2). 
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Fig. 1. — Transit photometry for HD 149026, with 40 second resolution (bins of 100 images). 
The top panel displays the raw light curve and the middle displays the light curve corrected 
for the detector ramp as described in section 3. At bottom are the residuals from the best 
fit light curve. 



- 7- 



0.0100 



D 

D 
X) 
'ifl 

CD 

X) 

c 0.0010 



o 



0.000 




10.0 



10 100 1000 
Bin Size 



10000 



o 

Q_ 

"O 
CD 
N 




0.0 0.2 0.4 0.6 0.8 1.0 

Frequency [min"^] 



1.2 



Fig. 2. — Left: Root-mean-square of binned residuals vs. bin size. The solid line is pro- 
portional to N~^/'^ and is normalized to match the value for bin size iV = 1. Right: Power 
spectrum estimate for the time-series residuals. The estimate has been divided by the power 
spectrum expected for randomly generated white noise of the same standard deviation and 
time sampling as the Spitzer data. The dashed line represents the median peak value of the 
simulated power spectrum estimates. 
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Because of the light curve modehng degeneracy between the parameters a and i?*, 
transit photometry alone cannot determine the quantity of interest, Rp. To break this 
degeneracy, one can either apply an external constraint on a (typically via Newton's version 
of Kepler's Third Law and a constraint on M^), on R^, or on some combination of both. 
Before applying any such constraints, we estimated the probability distributions for the 7 
light curve parameters by using the widely employed Markov Chain Monte Carlo (MCMC) 
technique (see e.g., Tegmark et al. 2004; Winn et al. 2007; Burke et al. 2007). The benefit 
of performing this analysis without any a priori assumptions on stellar quantities is that 
we can compare, on equal footing, the light curve constraints derived from our data with 
constraints derived from other photometric data. For this analysis, we adopted a conventional 
function as our goodness-of-fit statistic: 



where fmod{i) the calculated flux at the time of the i*'^ data point, fohsi}) is the flux 
measurement, and a is fixed to the value of the rms determined in the previous paragraph. 
Wc produced 10 chains of length 10^, with each chain starting from independent parameter 
points randomly chosen from a broad region (spanning approximately 5o") in parameter 
space. The beginning 25% of each chain was trimmed and the 10 chains were concatenated. 
We found that the Gelman- Rubin R statistic was < 1.01 for each parameter, which is an 
indication of convergence. 

Our MCMC analysis yields Rp/R^, = 0.05158 ± 0.00077 and impact parameter \b\ = 
|acosi/i?*| = 0.62 toM- Our result for Rp/R^ is larger than that of W08 {Rp/R^ = 
0.0491 laooos)' though the difference is only modestly significant. Our result for |6| is 
higher than that of W08 (|6| = 0.00 with 68% upper limit |6| = 0.36). We expect that 
our Rp/ Ri, result is more robust than the result from the optical data because, for weakly 
limb-darkened IR light curves, the radius ratio is measured almost directly from the observed 
fiux decrement. For optical light curves, however, the radius ratio is strongly covariant with 
the assumed limb- darkening coefficients. That is, an error in the assumed limb-darkening 
can translate into an error in the radius ratio. To investigate this point, we conducted a 
comparison with previously published {h + y)/2 HD 149026 light curves of Sato et al. (2005) 
and W08 (see §4.2 for further discussion of these data). For this sub-study we modeled the 
data in a manner similar to the above, but assumed a linear limb-darkening coefficient, which 
we allowed to vary freely. For the optical light curves, we determined a correlation coefficient 
between Rp/R^^ and the limb- darkening coefficient of r = 0.55, while for the same experi- 
ment with the Spitzer data wc found a much weaker correlation (r = —0.20). Furthermore, 
in the presence of strong limb-darkening, the radius ratio and impact parameter are also 
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correlated; the radius ratio can be traded off with the impact parameter to produce similar 
transit depths. We point this out because it suggests that the above mentioned discrepancies 
for b and Rp/ R^, are in fact correlated with each other. 

The results for these parameters and other important transit observables are reported 
in Table 1 (marked with a superscript 'a'). Noteworthy are the results for the mean stel- 
lar density, p^, and planet surface gravity, gp, which are model-independent determina- 
tions making use of information only from transit photometry and Doppler measurements 
JSeager fc Mallen-Ornela^ I2OO3I : ISouthworth et al.lboOTi : ISozzetti et aPboOTh . We also find 
a/R^ = 6.23 toll, which is consistent with the determination of W08 {a/R^ = 7.11 ^0:81). 
Note that the corresponding fractional uncertainty in our result for Ri,/a is fairly large 
(~ 10%). 



4. Stellar and Planetary Properties 

In the transit modeling literature, the parameter of interest, Rp, is usually determined 
via one of the following two methods. In the first, one obtains an externally determined value 
of R^, and then multiplies by the light curve results for Rp/R^,. In the second, one assumes 
a value for M^,, utilizes Newton's version of Kepler's Third Law to derive the semi-major 
axis, a, and then applies the light curve results for -R^/a (and -Rp/-R^)EI While the first 
method has the advantage that transit photometry determines Rp/R^, more precisely than 
Rp/a, the resulting Rp depends strongly on the assumed i?* {Rp oc i?^). This method also 
has the disadvantage of effectively disregarding any information gleaned from the light curve 
on R^/a. In the second method, the result for Rp depends only weakly on the assumed 
{Rp oc M^^), but, in our case, Rp/a is not constrained well enough to lead to a satisfactorily 
precise determination of Rp. For these reasons, we adopt a hybrid approach, imposing a 
radius constraint and, to make use of the -R*/a information, a mass constraint. Though the 
addition of this mass constraint represents an increased dependence on stellar models, we 
consider it a fairly benign dependency given how weakly the mass enters into the transit 
modeling (cx mI^^). 

In this section, we augment the Spitzer dataset with 10 previously published light curves. 
Together with independent constraints on the stellar properties described below, we fit for 
the Rp and other planet and stellar quantities. 



^Another possible route is to assume a stellar mass-radius relation ( Codv fc Sasselovll2002f l. We have 
chosen not to make such an assumption because of the uncertainty in the age and evolutionary state of HD 
149026. 
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4.1. Stellar Radius and Mass 



Using a collection of interferometric angular diameter measurements, iKervella et al. 



(120041 ) derived empirical relations for the angular diameters of dwarf stars as a function of 
Johnson magnitudes. We use their V, K relation for the angular diameter, (mas). 



log0 = 0.0755(\/-K) + 0.5170 - 



(2) 



Kervella et al.l (120041 ) find the root mean square residual from this best fit relation to be less 
than 1% for 20 stars ranging from spectral type AO to M2. We applied V = 8.15, as found 
in the Hipparcos Ca talog, and K = 6 .85, after transforming the 2MASS Kg magnitude to 



Johnson K following ICarpenterl (120011 ). After propagating the uncertainties in the photom- 
etry and the Kervella et al. (2004) best-fit parameters, we determined = 0.1755 ± 0.0021 
mas (see also Torres et al. 2008). The formal uncertainty in angular diameter is thus 1% 
and negligible compared to the uncertainty in parallax. After combining with the re-reduced 
Hipparcos parallax and uncertainty (vr = 12.59 ±0.70) of Ivan LeeuwenI (|2007| ). we determine 
-R^ = 1.50 ± 0.09 R0. For the stellar mass, we adopt the value = 1.30 ± 0.10 M0 from 



Sato et al.l (120051 ). who derived the value by matching stellar evolution tracks to spectroscopic 
properties. 



4.2. Light Curve Analysis Revisited 

We simultan eously fitted our Spitzer data togethe r with 3 hght curves published by 



Sato et al.l (120051 ) . 2 light curves by ICharbonneau et al.l (120061 ) and 5 light curves by W08. 



The 10 previous transit observations are discussed in detail in the references above, but we 
describe them briefly here. The Sato et al. 2005 and W08 observations were obtained with 0.8 
m automated photometric telescopes at the Fairborn Observatory. Fluxes were measured 
simultaneously through Stromgren b and y filters and averaged to create {b + y)/2 fluxes. 
The Charbonneau et al. (2006) observations were obtained with the Fred Lawrence Whipple 
Observatory 1.2 m telescope through the Sloan g and r filters. In analysis of the 10 light 
curves, W08 divided each raw light curve by a linear function of time that was fitted to the 
out-of-transit data. This step corrects for airmass effects and other systematic trends, but 
also has the effect of normalizing each light curve to have unit mean out-of-transit flux. We 
adopted these corrected data as well as their revised photometric errors, which were rescaled 
to account for the effects of noise correlation on ingress/egress timescales. Note that the 
composite of these 10 light curves, when binned to 30 second resolution, shows roughly the 
same scatter (~ 0.9 mmag) as the the Spitzer data binned to the same resolution. 

Next, we revisited the MCMC analysis of §3. We modeled the light curves as before. 
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using the small-planet transit routine of iMandel fc Agoll (120021 ). For t he q and r b and 



data, we assumed linear limb-darkening, with coefficients as tabulated by I Clare tl (120041 ) for 
a 6250 K, logg = 4.5, and [Fe/H]=0.3 star. For t he (b + y )/2 d ata, we assumed a linear 
limb-darkening coefficient of 0.712, the average of I Claret J2000h b and y limb-darkening 
coefficients, following W08. We employed 9 free parameters: Rp/R^,, i, cq, ci, C2, P, Tc, M^,, 
and R^, where the ramp correction coefficients, q, apply only to the Spitzer data. Note that 
we fit for only a single mid-transit time and required the transits to be spaced at integral 
multiples of P. W08 found no significant deviations from predicted transit times, so this is 
a reasonable assumption. We modified our goodness-of-ffi statistic as follows: 
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(3) 



with the second and third term reflecting the above determined stellar radius and mass with 
errors assumed to follow normal distributions. Note that the mass and radius constraints 
have entered into the in a simple additive form, which is strictly valid only if the con- 
st raints were determ ined entirely independent of each other. In fact, the mass determination 
of ISato et al.l (l2005l ) makes use of the parallax, which implies that the mass and radius deter- 
minations have some level of intrinsic covariance. To examine the impact of this covariance, 
we repeated our analysis with the following trial goodness-of-ffi statistic. 
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(4) 



where we experimented with values of the correlation coefficient, pmr, between -1 and 1. For 
— 0.7 < Pmr < 0.7, we found the best-ffi parameters and error bars to be negligibly affected 
by the correlation, and we found the results were significantly impacted only when \pmr\ > 
0.9. Since we expect the covariance between the stellar mass and radius determination to be 
more modest, we conclude that our results are not impacted by neglecting the covariance. 

We conducted the analysis as before; we produced 10 Monte Carlo chains of length 10^, 
cut the ffist 25% of each chain, and then combined the chains. For each parameter the 
Gelman- Rubin R statistic was well within 1% of unity. In table 1, we report best-ffi values 
and uncertainties for various parameters. We take the best-ffi value to be the median of the 
MCMC samples, and for the uncertainties, we report the interval that encloses the central 
68.3% of the MCMC samples. 



We determine a stellar radius of 1.497±0.069 R©, which is moderately refined compared 
to its prior distribution (1.50 ± 0.09 R©). This refinement indicates that the combination of 
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Spitzer and ground-based data is able provide some statistical influence on the parameter 
estimation through the observational constraint on i?^/a. We note that the external stellar 
radius constraint reinforces the high impact parameter solutions favored by the analysis of 
§3 {Spitzer data alone). This arises because the radius constraint favors relatively large 
stellar radii, which, for the given observed transit duration, can only be accommodated by 
non-equatorial impact parameters. The planet radius is determined to be 0.755 ± 0.040 Rj, 
with an uncertainty that is reduced versus previously published determinations (for example, 
0.71 ±0.05 Rj as determined by W 08). The reduction is partly due to the smaller uncertainty 
in the revised Hipparcos parallax (Ivan LeeuwenI 120071 ) that we have adopted, and partly due 
to the combination of the Spitzer and ground-based data. 
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Fig. 3. — Joint posterior probability distributions for Rp and directly observable quantities, 
as estimated by the MCMC analysis described in §4. The contours mark the 68% and 95% 
confindence regions. Left: Joint distribution for Rp and the "transit depth," [Rp/R^y. 
Middle: Joint distribution for Rp and the transit duration, defined as the interval from 1** 
to 4*'* contact. Right: Joint distribution for Rp and the ingress (egress) duration, defined as 
the interval from 1** to 2"'^ (3'"'^ to 4*^*) contact. Note that the vertical axes in each panel 
are scaled so that they encompass roughly the same fractional variation. 
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In Fig. [3|, we look deeper into the observational constraints on the key parameter, Rp. 
The most directly observable quantities from a transit light curve are the transit depth, 
the total transit duration, and the ingress or egress duration; a measurement of these three 
observables is sufficient (at least in the absence of limb-darkening) to determine the more 
physical parameters Rp/ R^,, Ri,/a, and cosi. By examining the joint posterior distributions 
for the three observables with the parameter Rp, one can gain insight into the current obser- 
vational limitations on the precision of Rp. While each of the panels in Fig. [3] demonstrate 
covariances, the third panel {Rp and ingress duration) shows particularly strong covariance. 
Thus, the major limiting factor in reducing the uncertainty in planetary radius appears to 
be the ability to resolve the ingress duration. Unfortunately, constraining this quantity with 
ground-based photometry is complicated by lini b- darkening and the effects of systematics 



and correlated errors (e.g.. |Pont &: Moutoul 120071 ) 



4.3. Influence of Star Spots 

As with limb-darkening, an inhomogeneous surface brightness due to spots would impact 
both the depth and shape of the transit light curve. If the transit chord intersects a star 
spot, a positive "bump" will be introduced into the transit light curve, while if the transit 
chord is along an unspotted area of an otherwise spotted star, the transit would appear 
deeper (see e.g. Knutson et al. 2008, Beaulieu et al. 2008). The existence of star spots can 
be investigated by long-term photometric observations of the star, monitoring for periodic 
flux variations. 

Previously published APT data has shown HD 149026 to be photometrically stable to 
0.0015 mag, the limit of precision of the APTs (Sato et al. 2005). Further (6 -|- y)/2 out- 
of-transit observations have been obtained with the APTs, extending the dataset to over 
3 years (Winn et al. 2008; G. Henry, private communication). With this additional APT 
data, kindly shared with us by G. Henry, we have searched for evidence of star spot-induced 
variability. We computed the periodogram for the time-series (550 total flux measurements) 
in fine steps of the period for periods between 0.5 and 100 days. Examination of all prominent 
peaks in the periodogram reveals no evidence for any significant periodicities, and allows us 
to place an upper hmit on the peak-to-peak amplitude of any sinusoid (in the period range 
0.5-100 days) of less than 0.001 mag. Any spots at or below this level will have negligible 
impact on the transit light curve, especially given that the 1 a uncertainty in the transit 
depth for HD 149026b is 3%. 
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4.4. Refined Ephemeris 

The precise transit timing from §3, along with the fact that the Spitzer hght curve 
extends the time base-hne of HD 149026b transit observations, enables a significant refine- 
ment in the transit ephemeris. For the previously published transit observations, we adopt 
the transit times and uncertainties listed in Table 3 of W08. We fit the timing data to the 
equation 

T,{E) = T,(0) + ExP (5) 

where Tc is the transit time, E is the transit epoch, and P is the orbital period. We determine 
P = 2.8758887 ± 0.0000035 and ^(0) = 2454327.37211 ± 0.00047, with xV^/'dof = 0.564, 
with N^of = 10. In Fig. HI we show the transit time residuals for all published transits. 
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Fig. 4. — Observed minus calculated mid-transit times for HD 149026b. The calculated 
transit times are derived from the ephemeris in eq. (2). The estimates for the first 11 transit 
times are drawn from Table 3 of W08. 
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Discussion 



We have presented and analyzed Spitzer 8 /im transit observations of the HD 149026 
system. By incorporating previously published data, and adopting constraints on the stellar 
mass and radius, we improve the determination of the planetary radius to Rp = 0.755 ± 
0.040 Rj. Our measurement reinforces previous findings of the intriguingly small radius of 
HD 149026b. To place this result in context, models in which HD 149026b (with a total 
mass of 114 ±2 ) is composed purely of H/He require a radius greater than 1.1 Rj (e.g.. 



Burrows et al.l 120071 ). 

The implications of the small measured radius on the i nterior structure o f HD 149026b 



have been modeled by a number of authors . Most works (ISato et al.l l2005l : iFortney et al. 



20061 : llkoma et al.l l2006l : iBurrows et al.l 120071 ) have assumed that all of the heavy elements 
reside within the planet's core, although it was often stressed that this may not necessarily be 
the case. For instance, recent models of Jupiter's structure indicate that the r najority of its 



heavy elements are mixed within the H/He envelope (Saumon & Guillot 2004). iBaraffe et al. 



(120081 ) recently computed evolution models of HD 149026b and other planets and showed 
that if these heavy elements are distributed within the envelope, rather than all in the core , 
less are needed to obtain the same model radius at a given age. However, llkoma et al.l (120061 ) 
also explored this effect, and noted that enhanced metallicity of the H/He envelope should 
also lead to higher atmospheric opacity, which will slow the contraction, necessitating more 
heavy elements. 

The choices made by the modelers have been diverse, and many different atmospheric 
boundary conditions, assumed heat capacities of the heavy elements, and equations of state 
(EOSs) for the heavy elements have been explored. Interior heavy element mass estimates 
have generally ranged f rom 50-90 for a planet radius of 0.725 ± 0.05i?j. On the low 
end, llkoma et al.l (120061 ) find that a 35M(^ core would be necessary, if the planet cooled and 
cont racted in isolat i on, th en was brought to 0.042 AU at the present time. On the high 
end, iBurrows et al.l (120071 ) find 110 M^, if the planet has an atmospheric opacity 10 times 
larger than solar composition atmosphere models. For all of these models, uncertainty in 
the measured radius is more significant than the uncertainty in the system age. 

A full exploration of new evolution models, including the potential contribution of 
TiO/VO opacity, which may be present in the planet's visible atmosphere (Fortney et al. 
2006, Fortney et al. 2008, Burrows et al. 2008), is beyond the scope of this paper. Given 
the previous modeling efforts, together with uncertainties in atmospheric metallicity and 
opacities and the distribution of heavy elements within the planet, the 50-90 heavy ele- 
ment mass range is still likely to be correct, even for our modestly larger measured radius. 
We note that current estimates of the heavy element abundance of Saturn (which is simi- 
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l ar in mass to HP 14902 6b) and Jupiter range from 13-28 Mq^ and 8-39 M^, respectively 
( ISaumon &: GuillotI l2004j ) . Uncertainty in the composition of giant planets is the rule, not 
the exception. 



As has been stressed recently by iBurrows et al.l (120071 ). and others, constraints for any 
particular planet will remain uncertain, but with a large sample size of transiting planets at 
various masses, radii, orbital distances, and stellar metalli city, trends will emerge which will 



shed l ight on the formatio n and structure of these planets (iGuillot et al.ll2006l : iFortney et al. 
200?! : burrows et al.l[2007l ). 



We are thankful to F. van Leeuwen for providing parallax data for HD 149026, and to G. 
Takeda for discussions regarding the spectroscopic determination of stellar properties. We 
are especially grateful to G. Henry for sharing many seasons of photometric data. We would 
also like to thank an anonymous referee for specific and helpful recommendations. This work 
is based on observations made with the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. 
Support for this work was provided by NASA through an award issued by JPL/Caltech. 
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Table 1. Estimates of the HD 149026 System Parameters 



Parameter 


Median 


15.9*'* Percentile 


84.1** Percentile 




0.05158 


-0.00077 


+0.00077 


Rp/Ri, 


0.05147 


-0.00077 


+0.00076 


i [deg]^ 


85?4 


-1?9 


+2?5 


i [deg] 


85?3 


-0?8 


+0?9 


a/R.^ 


6.23 


-0.63 


+0.71 


a/R* 


6.20 


-0.25 


+0.28 


p*[g cm~^]'^ 


0.51 


-0.13 


+0.21 


p*[g cm"^] 


0.547 


-0.064 


+0.076 


log^p[cgs]^'^ 


3.18 


-0.09 


+0.10 


logfi-pfcgs]^ 


3.203 


-0.048 


+0.049 


P [days] 


2.8758887 


-0.0000035 


+0.0000035 


Tc [HJD]- 


2454327.37213 


-0.00050 


+0.00050 


Rp[Rj] 


0.755 


-0.040 


+0.040 




1.497 


-0.069 


+0.069 



^Determined from analysis of Spitzer data alone 
^Using K = 43.3 ± 1.2 m s-^ from Sato et al. (2005) 



